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ABSTRACT 

We have made comparative studies of ion anisotropy and high-energy variability 
of solar energetic particle (SEP) events previously examined by the Solar, Heliospheric, 
and Interplanetary Environment (SHINE) Workshop campaign. We have found distinctly 
different characteristics of SEPs between two large “gradual” events having very similar 
solar progenitors (the 2002 April 21 and August 24 events). Since the scattering centers 
of SEPs are approximately frozen in the solar wind, we emphasize work in the solar-wind 
frame where SEPs tend to be isotropized, and small anisotropies are easier to detect. 
While in the August event no streaming reversal occurred, in the April event the field- 
aligned anisotropy of all heavy ions showed sign of streaming reversal. The difference in 
streaming reversal was consistent with the difference in the presence of the outer 
reflecting boundary. In the April event the magnetic mirror, which was located behind the 
interplanetary shock driven by the preceding coronal mass ejection (CME), could block 
the stream of SEPs, while in the August event SEPs escaped freely because of the 
absence of nearby boundary. The magnetic mirror was formed at the bottleneck of 
magnetic field lines draped around a flank of the preceding CME. In the previous SHINE 
event analysis the contrasting event durations and Fe/O ratios of the both events were 
explained as the interplay between shock geometry and seed population. Our new 
findings, however, indicate that event duration and time as well as spectral variation are 
also affected by the presence of a nearby reflecting boundary. 

Subject headings: acceleration of particles — shock waves — Sun: coronal mass ejections 
(CMEs) — Sun: particle emission — interplanetary medium 
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1. INTRODUCTION 

1. 1. Significance of Investigating High-Energy Variability of SEP Events 

One issue of concern to the National Space Weather Program (see 
http://www.nsf.gov/pubs/2007/nsf07010/nsfD7010.ispl is the investigation of high-energy 
variability of spectral and composition characteristics of gradual SEP events. We seek to 
know the origin of that variability, and the way it relates to outstanding questions in 
solar-terrestrial physics. 

Examination of particle spectral and composition characteristics in various SEP 
events is widely used to investigate the injection, acceleration and transport of SEPs (see 
Reames 1999). Observed spectral characteristics, however, strongly depend on the 
longitude of the point at which the observer’s magnetic flux tube connects to the CME- 
driven shock (e.g., Reames et al. 1996). There exists an east-west asymmetry of intensity- 
time profiles of SEPs that western SEP events reach peak intensities earlier than eastern 
events (e.g., Cane et al. 1988) because of the Archimedean spiral structure of the 
interplanetary magnetic field (IMF). In order to diminish the influence of solar longitude 
effects, it is preferable to make comparative studies on SEP events that have very similar 
solar progenitors. 

Often two gradual SEP events having very similar solar progenitors show similar 
characteristics at ion energies less than ~10 MeV nucleon' 1 , but at higher energies may 
exhibit extreme differences in their characteristics including abundance ratios, event size, 
spectral shape, GeV-ion content, and event duration (Tylka et al. 2005). Carrying out 
comparative studies on these differences can improve our understanding on the 
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generation and propagation of SEPs, and benefit the forecasting of the space radiation 
environment. 

1.2. The Solar, Heliospheric, and Interplanetary Environment (SHINE) Workshop 
Campaign 

The SHINE workshop (see 

http://cdaw.gsfc.nasa.gov/SHINE Campaign/index.htmH encouraged a concerted, 
focused effort to investigate a few carefully selected “campaign” SEP events. Two 
famous examples of SHINE events are the 2002 April 21 and August 24 events examined 
in Tylka et al. (2005, 2006), and Tylka and Lee (2006). The April and August events had 
their flare locations at S14W84 and S02W81, and CME speeds at 2400 kms' 1 and 1900 
kms' 1 , respectively. In the both events the size of associated flares was nearly same. The 
solar wind speed and the transit time of interplanetary (IP) shocks were also comparable. 
In addition, both events were accompanied by the metric and DH type III and type II 
radio emissions (Tylka et al. 2006). 

At ion energies between ~0.5 and -10 MeV nucleon' 1 the two events had nearly 
same event-averaged Fe/O ratio. At higher energies, however, the ratio differed by nearly 
two orders of magnitude: in the April event the ratio fell to only -10% of the nominal 
coronal value, but in the August event the ratio rose to -6 times the coronal value (see 
Figure 1 of Tylka et al. 2005; 2006). In addition, at <100 MeV the proton intensity-time 
profile in both events showed similar evolution patterns, but at higher energies the 
duration of proton intensities was distinctly different (see Figure 1 1 of Tylka et al. 2005). 
The full width of half maximum (fwhm) of proton intensities in the April event was -20 
hr in contrast with -2 hr in the August event. 
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Tylka et al. (2005, 2006) and Tylka and Lee (2006) attributed the difference 
between the two events to the interplay between two factors involved in the ion 
acceleration by CME-driven shocks (see Figure 2 of Tylka et al. 2005): the evolution in 
the geometry of CME-driven shocks, which generally begins as quasi-perpendicular near 
the Sun but evolves toward quasi-parallel as the shock moves outward; and a compound 
seed population, typically consisting of suprathermals from the corona or solar wind and 
from small impulsive flares. The quasi-parallel shock (Lee 1983) remains in contact with 
a given group of magnetic flux tubes for an extended period and has small injection 
energy requirement, resulting in a long-duration of accelerated particle events and a seed 
population similar to solar-wind suprathermals. In contrast, the quasi-perpendicular shock 
has a short contact period and a high injection energy requirement, preferentially 
accelerating seed particles from flares within a short period. 

1.3. Importance of Anisotropy Analysis ofSEPs 

So far the high-energy variability of SEP events has been examined mainly by 
using ion composition or energy spectral data. In addition to the ion composition or 
energy spectral measurement, however, the ion anisotropy analysis is an independent way 
to examine the high-energy variability of SEP events. In fact, SEP angular distributions 
are more direct means to study particle transportation in the interplanetary medium 
(Reames & Ng, 2002). 

Recently, Tan et al. (2007) carried out an anisotropy analysis of gradual SEP 
events in the 2-8 MeV nucleon' 1 range by using the Low-Energy Matrix Telescope 
(LEMT) data on the Wind spacecraft (von Rosenvinge et al. 1995). So far the 
Wind/LEMT sensor has provided the best resolution of ion angular distributions in the 
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MeV nucleon' 1 range (see Reames et al. 2001). The analysis began by introducing the 
concept of the “rest” frame, in which the phase space distribution function of ions is 
assumed to be isotropic (Gloeckler et al. 1984). The velocity of the rest frame relative to 
the spacecraft frame is the ion bulk flow velocity V F that can be calculated from ion 
sectored count rate data at given ion energy. Since in the solar wind frames the first-order 
anisotropy vectors A ls can be easily deduced from V F (Forman 1970; Tan et al., 1992a), 
we started from the V p analysis in order to examine anisotropic characteristics of SEPs. 

In two large events (the 1998 September 30 and 2001 September 24 events) among three 
SEP events examined in Tan et al. (2007) the flow reversal of heavy ions was observed in 
the spacecraft frame, while protons kept their flow direction continuously. 

A potential explanation why only heavy ions reverse their flow direction is that, in 
the given MeV nucleon" 1 range, softening spectra at the local IP shock may provide 
mainly accelerated protons, but fewer heavy ions (see Desai et al., 2003, 2004; Tylka et 
al. 2005, 2006). Consequently, heavy ions predominantly come from early acceleration 
near the Sun, and propagate across 1 AU. Beyond 1 AU, there is evidence (Bieber et al. 
2002; Reames & Ng 2002; Tan et al. 2007) indicating the possible existence of a nearby 
reflecting boundary of SEPs because of the transient structure of IMF driven by 
preceding CMEs. In fact, downstream of the IP shock a magnetic mirror can be formed in 
the bottleneck of magnetic field lines draped around a flank of preceding CME (Tan et al. 
1992b; Bieber et al. 2002). The magnetic bottleneck plays the role of outer reflecting 
boundary of SEPs. Forward streaming particles encountering the boundary could be 
reflected back to 1 AU to enhance the backward stream of particles. 

1.4 IMF Configuration in the 2001 September 24 Event 
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Since the presence of transient reflecting boundary of SEPs should be traceable 
from the observations of IMF and solar wind, we will reanalyze these observations in the 
2001 September 24 event previously examined in Tan et al. (2007) over a longer period 
prior to the occurrence of the primary CME2. The time profiles of the strength B of IMF 
and the speed V sw of solar wind from Wind observations are plotted on the top panel of 
Figure 1 , starting from ~1 week prior to the occurrence of CME2. The largest jump in 
both B and V sw occurred at time t obs (Shock2), indicating the arrival of the IP shock 
(Shock2) prior to the primary CME2, whose launch time was September 21 10:16 (UT) 
(linear fitting, see http://cdaw.gsfc.nasa.gov/CME_list) as denoted by t laun (CME2) in 
Figure 1. 

Beside the primary CME2, there was a preceding CME1. The Wind spacecraft 
observed the IP shock prior to CME1 at t obs {Shock\) = September 23 09:18 (UT), prior to 
the launch of CME2. The observed solar wind speed after t obs (Shockl) was V swl ~ 600 

kms' 1 . Assuming that the average CME speed between the Sun and 1 AU was between 
V swl and 3 V swl , the launch time t laun (CME 1) of the preceding CME1 should be between 
September 20 12:00 (UT) and September 22 10:00 (UT). During this time interval there 
were 16 CME events observed. Fortunately, only one CME event had its linear speed 
greater than 500 kms' 1 . That CME event with its linear speed of 659 kms' 1 should be 
identified as CME1, whose launch time t, aun (CME \ ) was September 21 08:48 (UT) 

(linear fitting). 

Since the magnetic bottleneck can be formed by the magnetic field lines draped 
around a flank of the preceding CME (Tan et al. 1992b; Bieber et al. 2002), the 
morphology of field lines behind Shockl should be essential to particle transport. We 
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hence show the field of view (FOV) for the preceding CME1 as observed by the large 
angle and spectrometric coronagraph (LASCO)/C2 telescope on board of the SOHO 
spacecraft at September 21 10:54 (UT) on the middle left panel of Figure 1, where CME1 
was in the southeast quadrant with its central position angle (CPA) of 135°. 

Furthermore, the event-associated flare was located at S19E63. Assuming an axial 
symmetric expansion of the preceding CME1, on the middle right panel of Figure 1 we 
schematically draw the envelop of the CME material in the solar wind (the 
“interplanetary coronal mass ejection (ICME)”, see Cane & Richardson, 2003) as 
projected on the ecliptic plane at t = t laun (CME2), at that time the primary CME2 was 
still near the Sun, while the radial separation A r between the leading edge of CME1 and 
the Wind spacecraft was small (-0.34 AU). From the cartoon in Figure 1 it is seen that 
the Wind spacecraft was located on “ open ” field lines. Because of the high speed of 
CME1, however, the field lines draped around the western flank of CME1 would be 
compressed in the region between Shock 1 and CME1, leading to the formation of 
magnetic bottleneck that plays the role of reflecting boundary of SEPs. 

1.5 Effect of Nearby Reflecting Boundary on Particle Transport 

An intuitive speculation is that in analogy to building a dam in a stream, where 
the water level and storage time in the reservoir increase, the reflecting boundary that 
blocks the flux tube would increase the peak intensity and duration of high-energy 
particles inside the tube, leading to a larger particle fluence that is the particle intensity 
integrated over the SEP event period. It is noticeable that observational evidences indeed 
support this speculation. For example, Roelof et al. (1992) noted that an inner 
heliospheric “reservoir” of SEPs could be formed behind a magnetic structure that is 
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created earlier by a superposition of plasma disturbances that inhibit the escape of SEPs. 
Sarris and Malandraki (2003) found that the electron event occurring within a converging 
IMF structure exhibits a remarkably longer decay phase in comparison with an event 
occurring within a diverging IMF structure. In addition, Reames et al. (1996) and Ng et 
al. (2003) calculated the decay time of ions trapped behind the shock and found that it is 
independent of ion energies. Also, Kocharov et al. (2005) simulated the intensity of high- 
energy protons in a closed loop of IMF and found that inside the loop there is a nearly 
perfect exponential decay of proton intensities with the decay time being significantly 
longer than that predicted by the usual diffusion-convection model. 

Our speculation is also supported by the observation in the 2001 September 24 
event. Since the flow reversal of heavy ions (Tan et al. 2007) and the presence of the 
preceding CME (the top panels in Figure 1) are observable, according to our speculation 
the September event should have a high peak intensity and long duration of high-energy 
protons. As shown on the bottom panel of Figure 1, the GOES-8 proton data confirms the 
speculation. 

1.6. Questions to be Addressed in This Work 

We wish to add the anisotropy analysis of SEPs to the examination of both 2002 
April and August SHINE campaign events. Since their main characteristics were already 
reported in Tylka et al. (2005, 2006) and Tylka and Lee (2006), the reader is referred to 
these publications. 

In this work the first question we address is whether, in the solar wind frame, the 
field-aligned first-order anisotropy of ions is different between the April and August 
events. Note that in this paper the ion anisotropy is determined relative to the solar wind 
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frame, in which the scattering centers are approximately frozen and the ions will become 
isotropic in the absence of other influences. This is the reference frame in which to 
observe small anisotropies and their reversals. If the answer to our first question is 
positive, our second question will be whether the difference can be attributed to the 
transient characteristic difference of ion reflecting boundaries, which should be 
recognizable from ICME, solar flare, IMF and solar wind observations. Finally, our third 
question is what other high-energy characteristics of SEPs would be affected by the 
difference of ion reflecting boundaries. 

Data from the Wind, IMP, and GOES spacecraft are used in this work. We first 
present the observed data on the ion first-order anisotropy, high-energy proton duration, 
and Fe/O ratio in both events. Then, we discuss the formation of a transient reflecting 
boundary and its implication on changing high-energy characteristics of SEPs. 

2. OBSERVATIONS 

2.1. Observed Data 

Beside the Wind and IMP 8 spacecraft data used in our previous work (Tan et al. 
2007), the high-energy proton data in the NO A A/GOES spacecraft (see 
http://spidr.ngdc.noaa.gOv/spidr/l are also used in the present work. Unlike data from the 
Wind and IMP 8 spacecraft, the GOES sensors provide no information on ion flow 
directions. 

Since the power-law index y of the ion phase-space distribution function is 
involved in the calculation of the first-order anisotropy of ions in the solar-wind frame, 
A ls , knowledge on ion energy spectra is necessary. While the energy spectrum of heavy 
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ions is obtained from Wind/LEMT data, the proton spectrum is deduced from IMP 8 
measurements. In the absence of IMP 8 data we have developed the technique to deduce 
the proton spectrum from NOAA/GOES observations. 

2.2. Intensity-Time Profiles of Ions 

Intensity-time profiles of ions in the April and August events are respectively 
shown in the top panels of Figures 2 and 3, where the dashed vertical lines indicate the 
launch time of the primary CME2 ( t laun (CME 2)) and the observation time of the IP shock 
(Shock2) that was driven by CME2 ( t obs (Shock 2)). In addition, in the April event from 
the near-Earth ICME table given in Cane & Richardson (2003) it is seen that a magnetic 
cloud (MC) appeared between April 20 00:00 (UT) and April 21 18:00 (UT) as denoted 
in Figure 2 by the shaded green region. 

From the omnidirectional intensity data of ions we first calculate the ion 
differential energy spectrum, from which the logarithmic differential intensity ( log( J m ) ) 

of ions given at the mean energy T m of selected energy channel is obtained. Upstream of 
IP shocks for both April and August events we observe T m ~3.6 and ~6.1 MeV nucleon 1 

for the Wind/LEMT low-energy (LE) (T = 2.5-5 MeV nucleon 1 ) and high-energy (HE) 

( T = 5-8 MeV nucleon' 1 ) ion channels, respectively. For them the lack of an intensity 
peak at the time of IP shock passage suggests that the effect of ion acceleration at 1 AU 
by local IP shocks was insignificant. In addition, in the April event additional intensity 
enhancements of ions were seen out of the MC region, in particular for protons. 

2.3 Bulk Flow Velocity of Ions 

The bulk flow velocity V F of ions (Tan et al. 2007) is a measurement of ion 
anisotropy in the spacecraft frame. On the 2 nd and 3 rd panels of Figures 2-3 we plot the 
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time profiles of its magnitude V F and longitudinal angle (p F , respectively. From the event 
onset the decrease of V F with time is generally seen for all ions. The time variation of <p F , 
however, is significantly different between protons and heavy ions. Since the ion angular 
distribution measured on Wind/LEMT is relative to the longitudinal angle <p B of IMF, 
instead of <j> F we plot the angular difference <p F - (j) B on the 3 rd panel. In the April event V F 
of both protons and heavy ions began along +B. Then heavy ions turned their V F to the 
-B, while protons kept their +B direction. In contrast, in the August event both protons 
and heavy ions kept their V F along -B direction. 

While the obvious difference of V F between the both events can be attributed to 
the presence and absence of a nearby reflecting boundaiy of SEPs in the April and 
August events, respectively (Tan et al. 2007), the details of V F reversal in the April event 
bothers us. It is seen that O and Fe ions at given ion velocity (i.e., energy per nucleon) 
reversed their V F at April 22 -06:00 (UT), which should be relevant to the flow reversal 
of heavy ions as observed in Tan et al. (2007). In contrast, both He ions at same and 
higher velocity showed a magnitude minimum and directional reversal of V F at -1/2 day 
earlier (April 21 -18:00 UT), indicating the possible effect of magnetic discontinuity on 
V F occurring at the MC boundary. In fact, during the 1998 May 2 MC event at the MC 
boundary a nearly zero field-aligned component of ion first-order anisotropies was also 
seen by Torsti et al. (2004) for 17-22 MeV protons (see their Figure 1), and by us for a 
few MeV per nucleon ions (data not shown here). To avoid complications near the MC 
boundary, in our further examination on the April event we will concentrate on later 
times when the Wind spacecraft was out of the MC. 

2.4. First-Order Anisotropy of Ions in the Solar Wind Frame 
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It is interesting to note that, as one kind of ICMEs, the MC should be also driven 
by the solar wind. Similarly, the effect of magnetic discontinuity at the MC boundary 
would be clearly seen in the solar wind frame. That raises another reason to examine the 
ion first-order anisotropy in the solar frame. Based on the bulk flow velocity V F 
measured in the spacecraft frame, we have developed the technique to deduce the ion 
first-order anisotropy A ls in the solar wind frame. According to equation (4) in Tan et al. 
(2007), we have 

A 1s =Y(V f -V sw )/v, (1) 

where v is the ion speed, y is the power-law index of the ion phase-space distribution 
function, and V sw is the solar wind velocity. From equation (1) it is seen that at V F » V sw 
A ls **yV f /v, while at V F « V sw A u is antisunward. 

For the April and August events, the magnitude A u and azimuthal angle <p ls of 
A ls are shown on the 4 th and 5 th panels of Figures 2 and 3, respectively. In both events 
during the “onset” phase we have V F > V sw and during the “plateau” phase we have 
V F <V SW , where the classification of event phases is according to Lee (2005). 
Consequently, during the onset phase A ls showed a forward streaming along B, while 
during the plateau phase A ls was along -V sw , toward the Sun. 

Since in the August event (see Figure 3) the onset of MeV per nucleon ions 
appeared significantly later than the launch time of the primary CME2 (i.e., 
tiaun(CME 2)), during the first few hours in the event we indeed detected the background 
particle stream left by previous SEP events. In addition, since in the August event the 
polarity of B was mainly sunward (see 5 th panel of Figure 3), the Y-axis on the bottom 
panel of Figure 3 is chosen to be -A llJIS . As a result, the positive Y-axis on the bottom 
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panels of Figures 2-3 always points antisunward. Thus during the onset phase in both 
events all ions showed outward streaming. The situation, however, was different during 
the plateau phase of corresponding ions. In the April event while protons had no field- 
aligned stream ( A 1jJIB ~ 0), all heavy ions showed sunward streaming. A nearly zero 
stream of ions was reached at the boundary of MC, which was shown in Figures 2, 4, and 
5 as the shaded green region. In contrast, in the August event both protons and heavy ions 
kept their forward streaming away from the Sun. 

The two-dimensional A ls vector measured by the Wind/LEMT sensor can be 
decomposed into two components A 1jJ1B and A Ij1B , parallel and perpendicular to the 
projected component of B on the ecliptic plane, respectively. Since during the upstream 
period between t laun (CME2) and ? oii (Shock2) in the April and August events the 
magnetic field B was nearly in on the ecliptic plane (in the GSE system the mean 
latitudinal angle of B was 10 ± 20° and -18 ± 14°, respectively), these components can 
approximately represent the projected components of the first-order anisotropy vector 
along the directions parallel and perpendicular to B , respectively. Note that the positive 
direction of A UJIB is given along B, and the positive direction of A U JJS shows a +90° 
(anticlockwise) angle from that of A 1iJIB . The time profiles of A 1j1B and A UJIS for the 
April and August events are shown on the 6 th and bottom panels of Figures 2 and 3, 
respectively. It is seen that both events had similar A U1£ profiles, which could be caused 
by north-south density gradients and/or perpendicular diffusion (Zhang et al. 2003). 

These effects are generally difficult to be estimated because of the unknown diffusion 
tensor and particle density gradient. Since we are concerned mainly with the ion 
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anisotropy originating from ion streaming, the emphasis of our examination is put on 
A 1jJI8 , the field-aligned component of A ls . 

From the bottom panel of Figure 2 it is seen that inside the MC the field-aligned 
antisunward anisotropy A IjJI8 of all ions decreased with time. Like in the 1998 May 2 MC 
event (Torsti et al., 2004) we observe a nearly zero A WI8 at the MC boundary. 
Furthermore, out of the MC the A 1jJI8 of all heavy ions changed to be sunward, while 
protons kept a nearly zero A Ij1I8 . Therefore, by comparing Figure 2 with Figure 3 we see 
that in the solar wind frame during the plateau phase of corresponding ions the main 
difference between the April and August events on particle streaming is that the 
streaming reversal of heavy ions was observed in the April event, but not in the August 
event. In Discussion section we will consider other characteristic variation of A UII8 
during different phases. 

2.5 IMF and Solar Wind Observations 

Similar to the analysis in the 200 1 September 24 event shown in section 1 .4, here 
we examine whether the IMF and solar wind data are favorable to the presence and 
absence of a nearby reflecting boundary of SEPs in the 2002 April 21 and August 24 
events, respectively. On the top panels of Figure 4 we plot the time profiles of B and V sw 
in the April event as measured by the Wind spacecraft during ~1 week prior to the launch 
of the primary CME2 (i.e., t lmm (CME2)). Before t laun (CME2) there were at least two IP 
shocks observed at r ofa (Shockl-l) and / oiJ (Shockl-2). The CME corresponding to the 
first shock, whose observation time was at April 17 1 1:02 (UT), would have less effect 
on the April event because of its distant location. We hence only consider the CME 1-2 


that was relevant to the Shock 1-2. 
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Also, similar to what we did for the 2001 September 24 event, here we assume 
that the average speed of preceding CME1-2 between the Sun and 1 AU was between 
V swl (~600 kms' 1 ) and 3 V^,. Thus the launch times of CME1-2 should be between April 
16 1 1 :00 (UT) and April 1 8 09:00 (UT). During the suggested time interval there was an 
obvious candidate event of CME1-2 with its onset time at April 17 07:50 (UT) (linear 
fitting). The candidate CME was a fast (1240 kms' 1 ) halo event associated with both a 
M2.6 X-ray flare at S14W34, near the central meridian, and a gradual SEP event. In 
addition, from the difference between t laun (CME 1 -2) and r oi>J (Shockl-2) the estimated 
average speed of CME 1-2 between the Sun and 1 AU was 856 kms' 1 , being very close to 
that in the 2001 September 24 event. 

As shown in the top panels of Figure 4, the onset of the primary CME2 occurred 
inside the MC. Therefore, during the onset phase the field lines passed the Wind 
spacecraft were probably or mostly closed. Upon exiting from the MC, however, the field 
lines were open completely. Because of the complexity caused by MC boundary crossing, 
for the 2002 April 24 event it is difficult to draw a cartoon of field line configuration (like 
we did in Figure 1 for the 2001 September 24 event). Nevertheless, the only real 
questions are whether there were any reasonable observed driver CMEs to associate with 
the interplanetary shocks. The CME details are not relevant here. We hence expect that 
between the preceding CME and the IP shock prior to it there would exist a magnetic 
mirror that played the role of the outer reflecting boundary of SEPs. 

Finally, we briefly mention IMF and solar wind observations in the August event, 
whose B and V sw data are shown on the bottom panels of Figure 4, where the gap in 
Wind data has been filled with the OMNI combined data (http://cdaweb.gsfc.nasa.govA . 
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There was no peak of B with B m s 15 nT occurred until ~5 days prior to t laun (CME 2) , 
indicating the absence of nearby reflecting boundary in the August event. Nevertheless, a 
distant reflecting boundary of SEPs may still exist as evidenced from the observation of 
the IP shock (Shock 1) prior to the primary CME (CME2) as shown on the bottom panel 
of Figure 4. 

3. DISCUSSIONS 

3. 1. Duration of High-Energy Proton Intensities 

Intensity-time profiles of high-energy protons measured by the GOES-8 
spacecraft are shown on the top panels of Figures 5 and 6 for the April and August 
events, respectively. From the panels we note the following facts. (1) The time to reach 
the intensity maximum in the August event was shorter than that in the April event. (2) 
After passing the maximum, the proton intensity in the August event showed a relatively 
fast decay with the decay rate being both time- and proton energy-dependent. In contrast, 
the proton intensity in the April event showed a relatively slow decay with an exponential 
decay rate being independent of proton energies, which is consistent with the calculation 
of Reames et al. (1996) and Ng et al. (2003) behind IP shocks and the simulation of 
Kocharov et al. (2005) in loop-like MCs. (3) In the April event an intensity enhancement 
of <15 MeV protons was seen outside the MC region, while >15 MeV protons showed a 
continuous exponential decay with time when the MC boundary was crossed. (4) Similar 
to Tylka et al. (2005), we observe that the duration of high-energy proton intensities in 
the August event was much shorter than that in the April event. Therefore, the observed 
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difference of high-energy proton durations between the both events is in support of our 
speculation shown in section 1.5. 

3.2. Characteristic Variations of Field-Aligned First-Order Anisotropy of Ions in the 
Solar Wind Frame 

In the region upstream of IP shocks the A ljJlfi data in the April and August events 
are plotted on the 2 nd panels of Figures 5 and 6, respectively. From Figure 5 it is seen that 
in the April event during the onset phase A ljllfi was approximately ion velocity- 
dependent, because of the almost overlapping of A liJIS plots for various ion species given 
at same ion velocity. Since the diffusion model (Parker 1963) predicts an ion velocity- 
dependent field-aligned anisotropy, we compare the model prediction with observations. 
In the radial (r) diffusion model with the mean free path A = A ^ , where A,, and f 
being constant, for an impulsive release of ions at time t = 0 and r = 0 the predicted 
field-aligned first-order anisotropy is (see Appendix C of Ng et al. 2003) 

A 1s ,b = 3r /[(2 - /3)vfj . (3) 

Taking into account of the Archimedean spiral structure of IMF we substitute r = 1.15 
AU into equation (3) in order to calculate the time profiles of A UB given at the mean 

energies T m ~ 3.6 and ~6.1 MeV nucleon' 1 for the Wind/LEMT LE and HE ion channels, 
respectively. Here (3 = 0.4 ± 0. 1 is estimated from the comparison of high-energy proton 
intensities predicted by the diffusion model with GOES-8 observations in the August 
event (also see Ng and Reames 1994). The A lsllB values predicted by the diffusion model 
are shown on the 2 nd panels of Figures 5-6 as the dotted lines. 

In view of the fact that no free parameter is introduced into equation (3), during 
the onset phase in the April event the consistency between the prediction and observation 
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should be acceptable. During the plateau phase, however, the prediction and observation 
are opposite in phase, indicating the streaming reversal of heavy ions in the solar wind 
frame. In addition, for protons we see that A 1jIIB ~ 0, implying a balance between the 
forward stream of protons freshly accelerated by the IP shock and the backward stream of 
reflected protons earlier accelerated near the Sun, or, perhaps, it only implies a uniform 
radial intensity distribution of protons throughout the flux tube. 

Moreover, during the plateau phase in the April event we observed that the A lsB 
value was different among different ion species given at same ion velocity. We are 
interested in exploring the nature of such difference. Thus we plot the 
A u , B (He)/A u , B (0) and A lsi[B (Fe)/A UHB (0) ratios on the 3 rd and 4 th panels of Figure 5, 
respectively. Note the scarcity and scattering of data points during April 21 14:00-22:00 
(UT) when the MC boundary was crossed. In order to find the mean value of the above 
ratios during the plateau phase we calculate their weighted average over the period of 
2002 April 22-23 (UT). The deduced A Ui[B (He)/A uw (0 ) = 1.01 ±0.04 and 
A UflB (Fe)/A lsW (0) =1.2 ± 0.1. Since He and O ions with nearly same rigidities had same 
Aijis value, and the A 1jJIB value of higher-rigidity Fe ions was greater than that of lower- 
rigidity O ions, we conclude that in the April event during the plateau phase A 1j IIB was 
ion rigidity-dependent. 

The situation is different in the August event. Because of the delayed launch of 
MeV per nucleon ions we cannot examine its onset phase in detail. During the plateau 
phase, however, the prediction of the diffusion model was consistent with the 
observation, either in polarity or magnitude of A 1jIIB . The A lsllB (He)l A unB (0) and 
A^ B (Fe)/A ls liB (0) ratios are also shown on the 3 rd and 4 th panels of Figure 6, 
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respectively, although no firm conclusion can be extracted because of poor quality of data 
resulting from the extremely small values of A llJls . 

3.3. Time Profiles of Fe/O Ratios 

We calculate the Fe/O ratio over different ion velocity (energy per nucleon) 
ranges. We first parameterize the energy spectrum of ions by using polynomial fitting, 
from which the ion intensity integrated over a given velocity window is deduced by 
numerical integration. From the integrated ion intensities we then calculate the Fe/O 
ratio. Our deduced Fe/O ratio at T= 2.5-5 MeV nucleon' 1 is shown on the bottom panels 
of Figures 5 and 6 for the April and August events, respectively. 

In both events starting from the event onset, the Fe/O ratio presented an 
exponential decline with time, although the detail of time variations was different. The 
Fe/O ratio in the April event had a minimum (Figure 5), as predicted by Ng et al. (2003), 
while the ratio showed a monotonic decrease in the August event (Figure 6). Why would 
the Fe/O ratio enhance again after passing its minimum in the April event? We believe 
that the enhancement of Fe/O ratios could be relevant to the observed 
A UfiB (Fe)/A u , B (0)> 1, which implies that along the magnetic field direction Fe ions 
have a greater sunward flow than O ions. Consequently, more Fe ions should appear in 
the backward stream, leading to an increase of Fe/O ratios during the plateau phase. 
Effectively, the sunward flowing ions reflected from the mirror retain the high Fe/O 
values seen early in the event. 

3.4. Comparison with Previous SHINE Event Analysis 

In this work we add the anisotropy data of ions to the analysis of two SHINE 
campaign events previously reported in Tylka et al. (2005, 2006) and Tylka and Lee 
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(2006). We are mainly concerned with the change of ion transport in the 2002 April 21 
event due to the apparent presence of a nearby reflecting boundary of SEPs. Through our 
analysis we have observed the streaming reversal of heavy ions in the solar wind frame, 
indicating the possible existence of a transient reflecting boundary of SEPs. Evidence 
gathered from heliospheric, IMF and solar wind observations indicates that the magnetic 
mirror located between the preceding CME and the IP shock prior to it forms the 
boundary. In the April event during the onset phase the field-aligned first-order 
anisotropy of ions in the solar wind frame is ion-velocity dependent, while during the 
plateau phase the reversed streaming of ions is ion rigidity-dependent. 

It should be admitted that the presence of a nearby reflecting boundary of SEPs 
would significantly affect characteristics of SEPs in the April event. For example, the 
reservoir effect caused by the boundary would increase the peak intensity and duration of 
high-energy particles, leading to a high particle fluence integrated over the entire SEP 
event in space. In addition, since any boundary would have finite cut-off rigidity, at 
sufficiently high velocity (energy per nucleon), Fe ions may freely escape from the 
boundary, while lower-rigidity O ions would still be reflected. As a result, the boundary 
could cause a decrease of Fe/O ratios at very high energies. However, this latter effect 
should be indistinguishable from similar spectral “knees” produced by rigidity-dependent 
trapping during acceleration. 

In the previous SHINE event analysis, the variation of event durations, spectral 
shapes and Fe/O ratios were attributed to the interplay between shock geometry and seed 
population (Tylka et al. 2005, 2006; Tylka and Lee. 2006). Particle reflection at a 
boundary beyond 1 AU is unlikely to influence shock acceleration close to the Sun and 
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will not affect the magnitude of the Fe/O ratio but will affect SEP event duration as well 
as the time and energy variation of Fe/O. Our work does not question the importance of 
seed populations and shock geometry in determining the Fe/O ratio. However, the new 
observation presented in this work suggests boundary reflection is important to the 
interpretation of SEP characteristics observed during the April event. 

4. SUMMARY 

We have made comparative studies of the 2002 April and August SEP events that 
had very similar solar progenitors but showed distinctly different high-energy 
characteristics of SEPs. Our main findings are as follows. 

(1) In the August event the field-aligned anisotropy in the solar wind frame showed no 
signal of reversal, while in the April event the streaming reversal of all heavy ions were 
observed during the plateau phase. 

(2) In the April event the field-aligned anisotropy of both protons and heavy ions in the 
solar wind frame was nearly zero at the boundary of the magnetic cloud. 

(3) In the April event a shock wave from preceding CME with a peak magnetic field 
strength of ~20 nT was observed within ~1.5 day before the launch of the primary CME. 
In the August event, however, there was no preceding CME within ~5 days before the 
primary CME. 

(4) In the April event the peak intensity and duration of high-energy protons were much 
greater than that in the August event. In addition, in the August event the decay time of 
high-energy protons was relatively short and proton energy-dependent. In contrast, in the 
April event the long decay time of high-energy proton intensities was independent of 
proton energies. 


23 


(5) In the April event the minimum of Fe/O ratios was consistent with a higher backward 
field-aligned anisotropy for Fe ions than for O ions during the plateau phase. 

We have been able to interpret these observations in terms of the presence of a 
nearby transient reflecting boundary that modified the properties of the 2002 April 21 
SEP event and an absence of such boundary in the 2002 August 24 SEP event. Particle 
reflection at a boundary beyond 1 AU is unlikely to influence shock acceleration close to 
the Sun and will not affect the magnitude of the Fe/O ratio but will affect SEP event 
duration as well as the time and energy variation of Fe/O. Our work does not question the 
importance of seed populations and shock geometry in determining the Fe/O ratio. 
However, the new observation presented in this work suggests boundary reflection is 
important to the interpretation of SEP characteristics observed during the April event. 


24 


ACKNOWLEDGEMENTS 

We appreciate the use of data provided by the NSSDC CDAWeb, NOAA, and the 
SOHO LASCO CME catalog, which is generated and maintained at the CD AW Data 
Center by NASA and the Catholic University of America in cooperation with the Naval 
Research Laboratory. We thank K. Ogilvie, A. Tylka, and X. Shao for many helpful 
discussions. Also, we thank the anonymous reviewer for his/her valuable comments. C. 

K. N. is supported under NASA grants LWS04-0000-0076 and SHP04-00 16-0024. 


25 


REFERENCES 

Bieber, J. W., et al. 2002, ApJ, 567, 622 

Cane, H. V., Reames, D. V., & von Rosenvinge, T. T. 1988, J. Geophys. Res., 93, 9555 
Cane, H. V., & Richardson, I. G. 2003, J. Geopys. Res., 108, A41 156 
Desai, M. I., et al. 2003, ApJ, 588, 1 149 

. 2004, ApJ, 61 1, 1 156 

Forman, M. A. 1970, Planet. Space Sci., 18, 25 
Gloeckler, G., et al. 1984, Geophys. Rev. Lett., 11, 603 
Gopalswamy, N., et al. 2004, J. Geopys. Res., 109, A12105 
Kocharov, L., et al. 2005, J. Geophys. Res., 1 10, A12S03 

. 2007, ApJ, 654, 1121 

Kahler, S. W. 2001 , J. Geophys. Res., 106, 20,947 
Lee, M. A., 1983, J. Geophys. Res. 88, 6109. 

Lee, M. A. 2005, ApJ, S158, 38 

Ng, G. K., & Reames, D. V. 1994, ApJ, 424, 1032 

Ng, C. K., Reames, D. V., & Tylka, A. J., 2003, ApJ, 591, 461 

Parker, E. N. 1963, Interplanetary Dynamical Processes, Interscience, New York 

Reames, D. V. 1999, Space Sci. Rev., 90, 413 

. 2000, in AIP Conf. Proc. 528, Acceleration and Transport of Energetic Particles 

Observed in the Heliosphere, ed. R. A. Mewaldt et al. (Melville: AIP), 79 
Reames, D. V., Barbier, L. M., & Ng, C. K. 1996, ApJ, 466, 473 
Reames, D. V., Ng, C. K., & Berdichevsky, D. 2001, ApJ, 550, 1064 
Reames, D. V., & Ng, C. K. 2002, ApJ, 577, L59 



Roelof, E. C., et al. 1992, Geophys. Res. Lett., 19, 1243 

Sarris, E. T., & Malandraki, O. E. 2003, Geophys. Res. Lett., 30(21) 2079 

Tan, L. C., et al. 1992a, J. Geophys. Res., 97, 1597 

. 1992b, J. Geophys. Res., 97, 179 

Tan, L. C., Reames, D. Y., & Ng, C. K. 2007, ApJ, 661, 1297 

Torsti, J. Rihonen, E. & Kocharov L. 2004, ApJ, 600, L83 

Tylka, A. J., et al. 2005, ApJ, 625, 474 

Tylka, A. J., et al. 2006, ApJ, SI 64, 536 

Tylka, A. J., & Lee, M. A. 2006, ApJ, 646, 1319 

von Rosenvinge, T. T., et al. 1995, Space Sci. Rev., 71, 155 

Zhang, M., Jokipii, J. R., & NcKibben, R. B. 2003, ApJ, 595, 493 


27 


FIGURE LEGEND 

Fig. 1 (upper) Time profiles of the strength B of IMF and the speed V sw of solar wind in 
the 2001 September 24 event, where t laun and t obs are the launch time of CMEs and 
the observation time of IP shocks, respectively; (middle) Field of view for the 
preceding CME1 as measured by the LASCO/C2 telescope (left) and the IMF 
configuration in the 2001 September 24 event (right, see text); (bottom) Intensity- 
time profiles of high energy protons as deduced from GOES-8 data in the 2001 
September 24 event. 

Fig. 2 Time profiles of the ion logarithmic differential intensity (log(7 m ) , where 

J m [ctn 2 s~ 1 sr“ 1 ( Me V /m)' 1 ]), the magnitude (V F ) of the bulk flow velocity V F relative 
to the spacecraft frame and the difference between the azimuthal angle (p F of V F and 
the azimuthal angle (p B of IMF ( (j) F - (f> B ), the magnitude ( A u ) and azimuthal angle 
(<p Als ) of the ion first-order anisotropy in the solar wind frame ( A ls ), the A ls 

components perpendicular and parallel to the projected component of B on the 
ecliptic plane, A u ^ and A bllfl for the 2002 April 21 event. 

Fig. 3 Same as Fig. 3, but for the 2002 August 24 event. 

Fig. 4 Time profiles of B and V sw for the 2002 April 21 event (upper) and August 24 
event (lower), respectively. 

Fig. 5 Time profiles of high-energy protons from GOES-8 observations, the field-aligned 
anisotropy A UB as compared with the prediction of the diffusion model, the 
^\s,uB(He)/A unB (0), A lsMB (Fe)l A uub (0), and J(Fe)/J(0 ) ratios for the 2002 April 21 


event. 


Fig. 6 Same as Fig. 5, but for the 2002 August 24 event. 
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